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Abstract 

The stationary ID Schrodinger equation with a polynomial potential V(q) of degree N is 
reduced to a system of exact quantization conditions of Bohr-Sommerfeld form. They arise from 
bilinear (Wronskian) functional relations pairing spectral determinants of (N + 2) generically 
distinct operators, all the transforms of one quantum Hamiltonian under a cyclic group of 
complex scalings. The determinants' zeros define (N + 2) semi-infinite chains of points in the 
complex spectral plane, and they encode the original quantum problem. Each chain can now 
be described by an exact quantization condition which constrains it in terms of its neighbors, 
resulting in closed equilibrium conditions for the global chain system; these are supplemented 
by the standard (Bohr-Sommerfeld) quantization conditions, which bind the infinite tail of each 
chain asymptotically. This reduced problem is then probed numerically for effective solvability 
upon test cases (mostly, symmetric quartic oscillators): we find that the iterative enforcement 
of all the quantization conditions generates discrete chain dynamics which appear to converge 
geometrically towards the correct eigenvalues/eigenfunctions. We conjecture that the exact 
quantization then acts by specifying reduced chain dynamics which can be stable (contractive) 
and thus determine the exact quantum data as their fixed point. (To date, this statement is 
verified only empirically and in a vicinity of purely quartic or sextic potentials V(q).) 



We study the 1- dimensional stationary Schrodinger equation with a real polynomial 
potential V(q) of degree N (> 2) on the real axis, taken in the rescaled form 



- + [V(q) + \}ip(q) = 0, V(q) = +q N + v 1 q N ~ 1 + v 2 q N ~ 2 + ■■■ + v N ^q. (1) 



We write the collection of coefficients as v — (vi, ■ ■ ■ , Wjv-i), absorbing any constant term 
of V into the spectral parameter A (here the sign-reverse of the usual energy). 

All standard analytical treatments of this problem lead to asymptotic results at best, 
i.e., to semiclassical or perturbative expansions which diverge factorially and are not 



though sophisticated techniques can sometimes convert those results to converging nu- 
merical outputs, those pathologies signal that such approaches may be mishandling the 
exact analytical structure of the problem ([!]). 

In this work we further develop an alternative and entirely exact method, which until 
now was specialized to quantizing the spectrum in the fully homogeneous case v = 
H, [Tj|. There, it produced self consistent, Bohr-Sommerfeld-like quantization formulae 
whose implementation (through iterative numerical schemes) apparently reconstructs the 
exact spectrum (convergent iterations are always observed, but this fact remain to be 
proved). The derivation of those results was quite indirect, via Borel transformations and 
quantum resurgence; this made calculations lengthy, hard to generalize, and pending on 
still incompletely established regularity properties of the Borel-transformed solutions (cf. 



The approach now to be described towards the same goal is much more direct, and 
immediately applicable to general (inhomogeneous) potentials. It only proceeds through 
a pair of auxiliary eigenvalue problems on a half-line, described in Sec.l; the spectral 
determinants of these problems exactly express the wave-function data at the endpoint 
through the identity (j38|), and thus they inherit a bilinear (Wronskian) functional relation, 
eq. (|43|) , through an analysis closely following Sibuya 0. Exact quantization conditions are 
then readily extracted for the corresponding Dirichlet /Neumann spectra, representing our 
analytical end result: eqs.fljSp (Sec. 2). As in the homogeneous case, these determine the 
unknown spectra only as fixed points of (explicit) iterations; the convergence of the latter 
is then a meaningful separate question, which here remains conjectural but is numerically 
probed in the remaining Sees. 3-4 on several test cases. First on even potentials, with 
the previous spectra directly giving the eigenvalues of eq.(Q) over the whole line (Sec. 3). 
Then, because eqs. (|38|) also supply the wave-function data at the endpoint of the half-line, 
the straight variation of this endpoint yields the solution of eq.(p]) itself; this allows us to 
also test the formalism on a (ground state) eigenfunction calculation (Sec. 4). 

1 Spectral preliminaries 

1.1 Polarized boundary conditions 

We first cast eq.(|l|) into two eigenvalue problems with asymmetric boundary conditions. 
Specifically, we restrict the problem to a half-line where V is confining: [0, +oo), we 



even Borel-summable in general 




(|; thm. 1.2.1 in 0). 
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keep the square-integrability at +00, while we put a Neumann, resp. Dirichlet, boundary 
condition at the finite endpoint q = 0. Both conditions then define self-adjoint operators 
H + , resp. H~ , which have purely discrete spectra we respectively denote {E 2 k} and 
{^fc+i} f° r k — 0,1,2,...; these admit an asymptotic expansion which is the Bohr- 
Sommerfeld formula reexpanded in descending fractional powers of the energy, 

53M?£~(A; + l/2), fc^ocinN, v = fi, fx - 1/N, (i - 2/N, ■ ■ ■ (2) 



def 1 , 1 7 def f pdq it 1/2 / 1 \ / /3 In 



and b^_ j/N (v) = a polynomial in the v r (f < j), e.g., b^_ 1/N 



2v x 
'nN 



in complete generality, the series depends on the sector (Neumann vs Dirichlet): b v = b^, 
but not until v = —3/2; whereas we will invoke those coefficients (and related ones, see 
Sec. 1.4) only in the leading range {y > — fi}, before this complication appears. (If V is 
an even polynomial, the two spectra coincide with the even vs odd parity components of 
the spectral problem on the whole line, hence their asymptotics coincide to all orders.) 

Concerning the exponent \x (the growth order), the forthcoming arguments often as- 
sume the property /i < 1 (i.e., N > 2) and always fj, ^ 1: the case N = 2 (singular) must 



be ruled out 12 



Then, up to sign, the eigenvalues read as the zeros of the Fredholm determinants 

A ± (A) = f J] (1 + A/S*); (4) 

, even 
k odd 

these infinite products converge to entire functions of order \x in the variable A, also entire 
in the parameters v ||. 

At fixed (X,v), let ip\(q) denote a recessive solution of eq.(P (i.e., a solution expo- 
nentially decreasing as q — > +00, which is unique up to normalization). Then Va(0) 
vanishes simultaneously with A - (A), (and likewise for an d A + (A)). However, the 

relations connecting them retain cumbersome factors unless each of A and i[) is suitably 
renormalized. The procedure for the determinants is the well-known zeta-regularization, 
but a separate rescaling of the wave functions also contributes in parallel to optimal 
simplification. 

1.2 The spectral (or functional) determinants 

By Mellin-transforming eq.(|2|), the spectral zeta functions (in analogy with |7J]) 

Z ± (s)= J2 E k s (and Z(s) = Z + (s) + Z"(s)) (Res>^) (5) 

. even 
k odd 



3 



are seen to extend meromorphically to lower Res values, and to be regular at integer 
points. An effective analytical continuation is brought by an Euler-Maclaurin summation 
formula which uses the asymptotic information (^): the regularized form 

zHs) = h )™ {Y.Ek s + l E K-l E rrz^ s+ 1 for ^odd n ( 6 ) 

K— >+oo <- Z Z , , — S -\- V) ' 

k<K {!/>!/()} v ' 

converges as soon as z/ can be taken below (Res). At the regular point s = 0, this gives 

zHo) = b j±\. (7) 

The point s = further serves to define spectral determinants, as zeta-regularized 
products: det H ± = f exp[— Z ± '(0)] and, by straightforward extension, 

D±(A) d J det( J ff ± + A) = expf-Z^O)] A ± (A). (8) 
This has the even more explicit Euler-Maclaurin representation 

log D ± (X) = lim { E \og(E k + X) + l \og(E K + A) (9) 

K ^ + °° k<K 2 

~\ E M&(log£*-i)} for fc,AT 



2 

{^>o} 



even 
odd • 



Crucially, we can and will form these zeta-regularized products for general sequences 
admitting asymptotics of the form @, including complex spectra [0 in which case the 
formulae are interpreted by analytical continuation from {A > 0, v = 0}. The {Ej,} are 
the zeros of the zeta-regularized product and they completely specify it (as opposed to 
a Hadamard infinite product, for which at least a normalizing factor has to be supplied 
independently); this normalization commutes with global spectral translations (unlike 
eq-@). 

1.3 Absolute WKB normalization for recessive solutions 

A recessive solution of eq.([l]) on R + admits an exact WKB representation |J (for A > 
- inf R + V) 

^x qo (q) = u(qy 1/2 e~ ^o u{q ' W where, for q -> +oo, u(q) ~ Tl(q) = (V(q) + A) 1/2 . 

(10) 

This normalization awkwardly depends on the base point go, so we seek a more intrinsic 
one with q = +oo: 

U<i)=<<ir l/2 e S « ° U{q ' )Aq ' ■ (11) 
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However, only (u — n)(q r ) is integrable at infinity whereas u(q) ~ n(g) ~ q N ^ 2 , so we add 
a further prescription, 

/ u(q')dq' ^ / (u - U)dq' + / (V(q') + A)- s dg'| s= _ 1/2 (12) 

J q Jq J q 

where the latter integral is now defined by analytical continuation from the domain 
{Res > 1/N, X > — inf V}. For the basic example V(q) = q , in the notations of 
eq.(|), 

X (g N + X)- S dq = r(g ~^ )r( ^ ) A- s+1/jV / + °V + A) 1/2 dg = -^-A". 

N r(s) jo 2sm7r/i 

(13) 

In general, however, J q +co (V + A) s develops a singularity at s = —1/2, which will 
explicitly affect many formulae hence requires further study. This singularity is immedi- 
ately seen to be independent of q (since J q (V + X)~ s dq' is finite for q, q" finite), and of 
A thanks to the regular right-hand side (at s = —1/2) of the functional relation 

f) r+00 r+00 

dxl {v + xrs = ~ s l ( v + x r {s+1) - ( 14 ) 

That singularity is therefore an intrinsic feature of the potential V alone. 

The desired analytical continuation of J q + °°(V + A)~ s can be performed leftwards from 
{Res > 1/iV} using the functional relation ([14]) in reverse, in complete analogy with a 
method exposed for spectral zeta functions in [13|. It follows therefrom that the only 
obstructions to regular continuation lie in the leading large-A behavior of J q (V + A)~ s , 
which decomposes into powers A~ s+P , the integration of which (as A~ s+p+1 /(— s + p + 1)) 
yields only simple-pole singularities. Consequently, a convenient regularization for the 
definition ( |I~2"D allowing such a pole at s = —1/2 is this finite part prescription, 



r+00 , , c r /-+oo 

/ (V + X)^^ lim / (v + xy 

Jo s=— 1/2 Jn 



R 



1/2 [Jq '* ' S + 1/2 J ^ 

with R = Res s= _i /2 f + °° (V + A)~ s = R(v), independent of (q, A). (16) 

Jq 

The new normalization (|TT| ) is then fully defined; it commutes with spatial translations. 



We can alternatively describe eq.(|Tl|) in the same framework |13| but using q as inte- 
gration variable instead of A. Thus, on the basis of the expansion 

(n<7) + A) 1/2 ~£/W forg^+oo, (<r=^, (17) 

$q°°(y + A) 1 / 2 gets specified as the indefinite integral consistent with the "symbolic 
integration" rules 

r+C °J°Aj ( +00 J°A„>M /-? CT+1 /(^+l) (<^-l) 

-log? (cr = -l) 



/■+OO , r+00 r+00 c r 

/ (^ + A) 1/2 ~£A,/ g'W, / g'W = { 

•/ o * J o Jo 
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Substitution into eq.flTTD then yields the asymptotic behavior 



^ A (g)~g-W 4+ ^exp{- £ /3 CT <f +1 / \<r + 1)} (q - +00) (19) 

{a>-l} 

(with no outer constant prefactor); this shows that the solution (0) exactly reproduces 
the subdominant solution defined in (ch.2, Sec. 6). 

1.4 Residues, residues 

A digression is needed here to better understand the residue R from eq.( [T6|) : by relating 
various asymptotic expansions we will find an identity, eq.(pO|), among several indepen- 
dently defined constants including R. 

a) We first consider the quantum partition function of the operator — d 2 /dg 2 + V(|<7|) 
on the whole line, namely 



hoc 



9(t) = £ e~ tEk . (20) 

fc=0 

This admits a t J, asymptotic expansion obtainable as in the case V(q) = \q\ N , but 
incorporating the lower-degree terms of V as perturbations: 

6(t) ~ E s -^~ v > HO, -1/ = -fi, -11 + + ^, ■ ■ ■ (21) 

(as with eq.(p|), this expansion is taken only to some finite order here). The leading 
nonclassical contribution to eq.(^TJ) again arises from the \q\ N term for which it is known to 



be 0(i M ); therefore, below this order, eq. fl2T|) also applies to (and is more easily computed 
from) the classical partition function, i.e. 



e d (t) = I ^ e -(p>+v(M)t = -!=/'" e- V{q)t dq. (22) 
^ 2 27r V7rt Jo 



-00 



We now Mellin-transform both partition functions ( P0f) and (|22"D term by term, as 

— / 6(t)e- x H s - 1 dt = J2(E k + \y s = Z(s,\), (23) 

T(S) V0 t^n 



k=0 



/ + °°g cl (t)e-^t- 1 dt= r ^ //-^ f +CO (V( q ) + X)-^dq^Z cl (s,X). (24) 

r(s) 7o r(s)v7r Jo 

On the quantum side (eq.(|23|)), this zeta function Z(s, A) |TB| relates to eq.(^) through 
Z(s,0) = Z(s), and to eq.® through 



exp[-dZ(s, \)/ds} s=0 = D(X) = D + (\)D-(X); (25) 
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whereas eg. p^) , drawn from the rightmost eq.(|22]), creates a classical picture in full 
analogy; in particular, the counterpart of eq.(^) evaluates using 

[-dZ cl (s, X)/ds} s=0 = 2 / (V(q) + X) 1/2 dq (26) 

Jo 



(in the sense of eq. ([15])) . 

(The parallelism between the normalization procedures @ for the determinants, resp. 
( [372] ) for the recessive solutions, is now fully clear: they involve the same analytical con- 
tinuations to s = of Z(s, A), resp. Z c \(s, A).) 

The small-t expansion ([H]) itself Mellin-transforms to a large-A expansion 

Z(s^)^J2~c-u^f-X- s+u (A-+oo) (27) 



also valid strictly above 0(A s M ) for the corresponding classical function (|24|) ; and the 
analogous treatment for the derivatives at s = [T3] likewise yields 



logD(X) ~J2c- v [-T(-v)\ u } with [-r(-0)A°] d =logA (28) 

also valid strictly above 0(A~ M ) for the analogous classical quantity ([261). 

Below, we will specially consider the coefficient c ; by the preceding arguments it has a 
wholly classical origin, and it is: the coefficient of t° in 6(t) or 6 c \(t); the coefficient of log A 
in logD(A) or 2j +oo (V(q) + A) x / 2 dg ; the coefficient of A" 1 ("residue") in J +co (V(q) + 
X)~^ 2 dq (by setting s = 1 in eqs.([M|) and (0) for Z cl ). 

b) The zeros of the above determinant -D(A) being {— E^}, their asymptotic behavior 
is given by the Bohr-Sommerfeld formula @. The following specifically holds before the 
latter begins to depend on the parity of k, which may occur for v < —3/2. Then, the 
two large-energy expansions (|2|), (|28| ) relating to the same function have to match || ||: 
their exponents {^} must coincide, as is already manifest, and their coefficients must be 
related, actually as 

b v = c- v /T(l-v) (i/ > -3/2). (29) 

This result at once describes (the leading part of) the expansion (|2|) and shows that it is 
of classical origin up to 0(A~ M ) (excluded). We are then going to focus on the quantity 
bo which contributes a constant shift term to the semiclassical series (0). 

c) In the purely classical expansion (]T7|) of (V(q) + A) 1 / 2 for q — > +oo, we select the 
coefficient of namely (3-\ (cf. [El, ch.2). By simple power counting, this coefficient 
is independent of A (for N > 2) and vanishes for iV odd. For iV even, it is the residue 
at q = oo of the complex-analytic function (V(q) + A) 1 / 2 ~ f q( N / 2 ^ (single- valued) near 
q = oo. We also denote it (3-\{v) to stress that it is an invariant of V. 

Our present goal is to obtain the following identifications, with R as in eq. (]T6|) : 

Z(0) ife { = } d ( = } -2R { = ] ~P-i(v) (= for N odd). (30) 
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Proof: (i) immediately follows from eqs.(|,|7]); (ii), from eq.(|29|) specialized at v = 0; 
(Hi), from comparing the residues at s = of the two integrals in eq . fl2~4|) ; finally, we 
prove (iv) by brute force. 

On the one hand, we extract R from eq.(|T6D but with q = and integration variable 
rescaled by A -1 ^, as 



(V+X)- s = \- s+1 / N / (a^+1) 



14 



4?„N-1 



l-N 

H h v N ^iA n x 



X N + 1 



ds. (31) 



The residue only affects the terms with overall weight A s 1 ^ 2 ; when the power in brackets 
is expanded, the corresponding coefficient comes as 



} N-1 



k>0 



{rj} 



(x N +l)- s - k x {N ~ 1)ri+ - +rN - 1 dx, (32) 



the inner summation being over (N— l)-uples {rj} subject to J2j=i r j — k and J2j=i j 
1 + N/2; hence the last written integral reduces to 



l/2)r(fc- 1/2) 
NT(s + k) 



(33) 



which acquires the residue 1/N at s = —1/2 for any k. 

On the other hand, we seek {3-\ from the large-g expansion 



V(q) 



1/2 



Q 



N/2 



VN-l 



1/2 



q q> 

when the power in brackets is expanded, the coefficient of q~ x comes out as 



•lwl ,1 . ^^V? 



} N-1 



E( 5 )(5+i)-( 5 -*+i)z : ri ,.. r ^, 

fa} 



(34) 



(35) 



fc>0 



with exactly the same summation range as before, hence this is manifestly N times the 
residue of eq . (|32|) at s = — 1/2. QED. 

Finally, again by simple power counting arguments, we remark that these classical 
invariants are not only absent for odd N, but also vanish for another broad class of 
potentials: all purely even polynomials of degree N multiple of 4 (as well as for all N in 
the special homogeneous case: /3_i(0) = 0.) Barring N = 2, nontrivial residues then first 
occur for non-even quartic potentials: 



0-i 00 



V 3 ViV 2 



4 + 16 



(V(q) = g 4 4- Viq z + v 2 q 2 + v 3 q), 



(36) 



and, within even potentials, for sextic ones: 



(V{q) = q 6 + v 2 q 4 + v 4 q 2 ). 



(37) 
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1.5 Basic identities 



Under the above notations, very simple identities connect the spectral determinants and 
the absolute-normalized solution: 

D"(A)=Va(0), D+(A) = -^(0). (38) 

(The proof is an adaptation of the arguments in ||, Apps. A and D). 

Next, following we continue eq.([TJ) in the complex g-plane down to the rotated 
half-line lying in the adjacent Stokes direction, namely [0, e -1 ^ 2 oo) where 

if = f — - (spectral symmetry angle). (39) 



By simple complex scaling upon q, 

^^^ X {j*' 2 q-M 1] ) (40) 
provides another solution (to eq.(|I|)) now recessive in the Stokes direction e _1</3 / 2 , where 

^1] def (eW 2 ^ e %/2 . . . ; e (7V-l)W2 ^ (41) 

expresses an action of the discrete rotation group of order (iV + 2) on the coefficients; 
equivalently it acts upon the potential V, mapping it to VW, then V®, ■ ■ ■ (now complex 
potentials). The order of the effective symmetry group is 

L = N + 2 generically, L = N/2 + 1 for an even polynomial V. (42) 

Now the Wronskian of the two solutions ip, ip^ of eq.([|), a constant, can be evaluated 
explicitly from their respective asymptotic forms (19), both valid as q — > +oo 0, and 
also expressed at q = by means of the respective identities fl3~8j ) for the potentials V and 
VM. Matching the two calculations then yields the fundamental bilinear identity 

e +i ^^ D + (e~^ X,v^)D-(X,v) - e-^ 4 D + (X,v)D-(e-^ X,^) = 2ie^- 1 ^/ 2 . (43) 

We stress that our approach will bypass any other matching of solutions, like those 
required in connection problems between nonadjacent Stokes directions, which yield non- 
trivial Stokes multipliers... 

Remarks: a) in the harmonic case (V(q) = q 2 , N = 2, ip = n), many of our argu- 
ments become invalid but the end result ( [i3|) remains, except that the residue is now 
A-dependent, /3_i = A/2 from eq.(]T7D (the identity ( f43|) now verifies by applying the re- 
flection formula for r(z) to D+(X) = 2^ 2~ x / 2 / T( 1 ± A ), D~{\) = ^ 2- A / 2 /r(^)); b) 
identities similar to eq.(El) have also surfaced in quantum integrable theories, as "quan- 



tum Wronskian conditions" 15 
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2 Exact quantization conditions 



For a homogeneous potential, we currently believe that the functional relation (43) and the 
asymptotic law ($) (imposed to some o(l) accuracy) together suffice to specify the whole 
spectrum exactly: indeed, we empirically recovered the spectrum as the fixed point of an 
(apparently) contractive mapping built using just that input f|, H OH- Thus, for zeta- 



regularized products built over sequences with prescribed Bohr-Sommerfeld asymptotics, 
the single Wronskian identity ( fi3|) can become so coercitive as to fix both its arguments 
completely. We now make the same guiding idea work for general polynomial potentials. 



2.1 The analytical result 

Guided by the homogeneous case ||, we strive to turn eq. (f43"D into equations where at 
least the Neumann (+) and Dirichlet (— ) spectra appear decoupled. To keep, say, only 
the former (the other admits a mirror-image treatment), we take eq. fl43|) and its partner 
with t/t -1 ! (written invoking the homogeneity property /3_i(t?W) = (?;)), then we 

set A = — E 2 k and eliminate D~(\, v) from the resulting pair, to find 

D + (e~ iv A, v [1] ) / D + (e +iv A, v^ 1] )\x=~E 2k = - e -'^ 2 + i ^-^ . (44) 

Whereas the procedure for V(q) = (^closed upon itself at once, this general one invokes 
the complex-rotated potentials and, step by step, all the L partner potentials V™ 

to reach closure. The spectra {E^l} (mostly complex) then altogether make up the 
independent unknowns (they are not independent as analytical continuations of each 
other in complex -y-space, but we are unable to get explicit relations expressing this). 

Next, absolute phases are specified in eq.(|3]) (and its cyclic permutations on the v$) 
by reference to the homogeneous case ||. This operation is essential to create quantization 
relations in the Bohr-Sommerfeld form (governed by an explicit quantum number k = 
0, 1, 2, . . .), but now exact. Our final result, in a synthetic notation, is an uncoupled pair 
of systems (using either (+, k even) for Neumann, or (— , k odd) for Dirichlet), each built 
of L exact quantization conditions: one per potential (with £ integer mod L), 

£g(if ] ) = 7r[fc + ~ ± *~ 2 . 1 + (-i)VMiO fo r k = lit: 



2 2(N + 2) 



where s|(E) d = -i[\ogD ± (- E,v [e+11 ) - \ogD ± (- e + '^ E,v [ "- ll )\, (45) 

the branches of logD ± (A) being taken by continuity from (A = 0, v — 0). 

The harmonic case V(q) = q 2 (with N = L = 2, ip = it) holds again with a degenerate 
structure, T^{E) = but = -E/2. 

2.2 Discussion 

All following considerations are meant in either spectral sector independently (Neumann, 
resp. Dirichlet), with all quantum numbers accordingly kept even, resp. odd. 
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Eqs.fl45|) form a system of constraints tying each E k at fixed £ to the two other spectra 
{Ej {Ef +1 ^} (whose zeta-regularized products build the determinants defining £±). 
The system of all such points in interaction is better displayed, in proper relative positions, 
if each spectrum is suitably rotated (we then call it a 'chain', as a more general name): 
for £ — 0, . . . , L — 1 (mod L), the £ th such chain is the set {e l£ip }, and each of its points 
is under the influence of every point in the two adjacent ((^± l) th ) chains — through the 
logarithm of their complex difference which enters eq.(^) via the formula (|9]). (In the 
generic case L = N + 2, adjacency is correctly shown on a double covering of the circle of 
asymptotic directions; the circle itself suffices only in the purely even case L = N/2 + 1.) 

The £ th equation is now a complex one precisely when its unknowns E$ are themselves 
complex, i.e., for £ ^ or L/2: eq. (flop thus remains a formally 'complete' system of 
mutual constraints for the unknowns E^. As in the homogeneous case we then surmise 
that eqs. (pE5|) are not only exact, but also genuinely complete, quantization conditions; i.e., 
they have the capacity to determine all their unknowns provided the asymptotic condition 
(0) is also enforced upon each spectrum {E^} separately (using the rotated coefficients 
v^). If moreover this resolution can be performed in any constructive way, then we may 
argue that the analytical formulae ( pP)| ) themselves are "giving" the solution of the stated 
spectral problem (Dirichlet or Neumann on the half-line). 

We will now report numerical experiments, mostly upon quartic potentials, which em- 
pirically confirm our conjecture in some regions of parameter space, by achieving effective 
computations of the spectra out of eqs. ([45|) (plus eqs.@). 



3 Numerical tests of exact spectrum quantization 

At present we can further resolve eqs.(^) only numerically. We will essentially seek to 
solve them by successive approximations for each level E$ in turn (by Newton's root- 
searching method), but stopping at some finite fc max , beyond which all eigenvalues can 
assume semiclassical values instead, once for all. The height of the cutoff k max , together 
with the depth in use for the semiclassical expansion (0), control the final accuracy. As for 
the zeta-regularized products involved, we use formula @ (again adding higher expansion 
terms to improve the A; max — > oo convergence). 

The preceding formalism immediately provides the quantization of levels for an even 
polynomial potential on the real line: parity symmetry splits the spectrum into even and 
odd sectors, which exactly correspond to the Neumann and Dirichlet problem on the 
half-line, respectively. As a side effect, parity also halves the order of symmetry L (cf. 
eq-©). 



3.1 Even quartic oscillators 

We now test the exact framework upon the quantization of levels for even potentials 
V(q) = 9 4 + v 2<f on the real line (the order of symmetry being L = 3). 
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To the initially real v 2 giving the real spectrum {Ek} are then associated: a complex 
spectrum {E' k } for the coupling constant jt>2, and {E^} for j 2 t> 2 (where j = e 27rl//3 ). The 
three chains {E k }, ]{E' k }, ] 2 {E k } are shown on fig.l for selected values of the harmonic 
coupling constant v 2 . They altogether constitute the reduced dynamical unknowns; by 
reality symmetry, there are only two independent chains (one real, and one complex). 

Then, in order at once to reach an equilibrium point for all constraints and to gain 
evidence for its uniqueness, we iteratively apply an elementary step: to recompute each 
chain in turn as (numerical) solution of its eq.(^) in terms of the adjacent chains at 
their current locations. Through an appropriate succession of such steps, we then try 
to build up a global process yielding a contractive iteration, as in the homogeneous case 
(now many more fine details can vary, and we tried but a few combinations). For us, 
a numerical validation of the exact quantization formulae results if one such iteration is 
found which displays geometric convergence to the correct spectra, i.e., a contraction ratio 
can be estimated with a minimum of numerical stability, and the limiting chains (checked 
upon their lowest five points, say) agree with independent eigenvalue calculations (e.g., 
matrix diagonalization). 

In order to define the mandatory asymptotic behaviors of all chains, we also need 
reference semiclassical chains for all three rotations of v 2 . These we take as the (numerical) 
solutions of eq.(^) pushed to 6 terms (=0(E~ 1 ^ 2 ), for better k — > +00 convergence). The 
same chains also conveniently serve as initial iteration data. 

Such an iterative approach was found for the homogeneous case (v 2 = 0) and it yielded 
strong contraction ratios (< 0.4) P, |J. However, this case also enjoyed full ternary 
rotation symmetry, and a qualitative change occurs now that this symmetry has to be 
relaxed. If we try to deform those earlier v 2 = iterations into ^-dependent schemes, 
the latter must apply synchronous updating to keep with the full ternary symmetry at 
v 2 = (i.e., the chains are recomputed individually but get actually updated all at once, 
at each completion of a full £-cycle). Unfortunately, such iteration schemes also appear to 
excite new (symmetry-breaking) chain fluctuation modes possessing much lower, almost 
marginal, stability (contraction ratios ~ ±0.9, down to v 2 = 0). As a consequence, 
their convergence discontinuously drops in the immediate vicinity of v 2 = 0, down to 
unacceptably low levels. On the other hand, we empirically found that iteration schemes 
of a different type, which immediately update the chains one by one, can remain fairly 
stable in a finite neighborhood of v 2 = (even though they are not as good as the earlier 
ones at the value v 2 = 0, where now the ternary symmetry is spontaneously broken). 

The simplest such scheme ("A") consists of recomputing the real and the complex 
chain alternatively. For v 2 near 0, this already appears to converge nicely and it indeed 
yields a first validation of the exact formalism. Still, its behavior deteriorates as v 2 
grows: between v 2 ~ +2 to +3, the iterations start to converge more erratically, then 
our root-searching algorithm (the Newton method) goes unstable, and shortly afterwards 
the contraction ratios of the iteration (plotted on fig.2, left) shoot up to almost unity 
(for v 2 < 2 they were estimated both from the iteration sequence and from diagonalizing 
the linearized dynamics near the fixed point, and for v 2 > 2 only the latter way). In 
particular, scheme "A" cannot be carried towards the harmonic limit. By contrast, it 



12 



behaves much better in the negative direction (double- well region), showing no sign of 
degradation over our test range, i.e., down to t>2 ~ —10 ! 

The positive-^ breakdown is possibly explained by fig. 1 (right): when t> 2 grows, the 
complex chain and its conjugate become almost degenerate towards low quantum numbers 
(a sort of complex tunneling effect), causing larger logarithmic pair interactions and in- 
stabilities (with the linearized-dynamics matrix entries growing like [mutual distance]" 1 ). 
An obvious remedy is then to remove the interaction between any two complex-conjugate 
chain points from the dynamics, by enforcing this symmetry as an a priori constraint 
(immediate updating also has to be kept); the so modified scheme ("B") exhibits a much 
more uniform linear stability indeed (fig. 2, right). This scheme converges and allows to 
validate the exact formalism for larger values of t>2, up to t>2 ~ +5. Beyond, we again 
become unable to get any convergence; so, global instabilities must still be creeping in, 
albeit more slowly. Again, the root-searching algorithm (the Newton method) diverges 
first, but now the linear contraction ratios themselves hardly grow at all (as obtained 
by diagonalization). Hence scheme "B" seems not to break down like "A" but rather to 
reach certain practical limits, beyond which only more elaborate implementations might 
resolve the case. (Perhaps, e.g., a more robust root-searching method would work, or the 
naive branch prescriptions implied in eq.fl45|) ought to be revised at such a distance from 
the homogeneous case.) 

The preceding break-points referred to iterations in the even-parity sector. As fig. 2 
shows, the odd sector tends to behave more stably. A conceivably better idea might thus 
be to confine iterations to the odd spectrum, then to solve for the even determinants from 
the coupling relations d43f) instead. 

In conclusion, the preceding tests favor the validity and effectiveness of the exact 
quantization mechanism for even quartic potentials over a sizable range of the harmonic 
coupling v 2 about 0, but are currently inconclusive for t> 2 > +5. 

In addition, we lately extended the scheme to sextic even potentials, with qualitatively 
similar results for the (still few) cases tested. 



3.2 Extension to non-even potentials 

In [|T2] we showed that the exact quantization formalism was fully valid for homogeneous 
potentials of any odd degree N, and specially for N = 1 which has Airy functions both 
as solutions and as spectral determinants. The paradox in the very regular behavior of 
this case is that the underlying potential V(q) = \q\ is not even once differentiable at 
q = 0. This example paved the way to the present generalization, establishing that the 
'even/odd' decomposition needed in the homogeneous case was just a Neumann/Dirichlet 
splitting, any parity properties of the polynomial V being irrelevant: an even potential 
over the whole line is always present as V"(|g|) (once its continuous differentiability is 
recognized as immaterial). 

We thus also tested the exact quantization conditions (B5j) upon a few non-even quartic 
polynomials (now for the Dirichlet /Neumann spectra). The exact quantization conditions 
then involve L = 6 distinct chains: £ = 0, 3 are real, (1,5), (2,4) being doublets of 
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complex conjugates. We found one iteration sequence ("C") to converge better than 
others (for no clear reason): I = {0,2,3,1} (then cyclically continued, with immediate 
updating of each chain and straight enforcement of 1-5 and 2-4 symmetry). The overall 
results were then comparable to those above for even polynomials (though occasionally 
achieving a lower numerical accuracy). 

This extension of the exact quantization formalism to non-even polynomials crucially 
opens the way to another application (and validation) now to be described. 

4 Exact wave- function analysis 

4.1 Analytical reconstruction of the wave function 

If no parity property is imposed upon the polynomials V(q) with respect to the endpoint 
q = 0, then this endpoint can be taken to an arbitrary value a (on the real line, at least), 
thus restoring translational invariance effectively. Equivalently, is kept as endpoint 
but the potential gets shifted to V a (q) — V(q + a) — V(a) on the half-line (becoming 
V(\q\ + a) — V(a) on the whole line). 

So, we now use eq. (|38D at q = a instead of and from right to left, to state that the 
absolute-normalized solution of eq.(|l]) at q — a, resp. its first derivative, are specified as 

Ma) = D-(V(a) + A), ^(«) = -D+(V(a) + A), (46) 

in terms of the spectral determinants for the modified potential V a (q) on [0, +oo), de- 
noted Df. But these in turn are canonically specified as the zeta- regularized products 
over their own chains of zeros; now the latter, as the eigenvalues of (V a )^(g), are given 
by exact quantization conditions of the form (fi5|), hopefully through some convergent 
iteration scheme as above. Thereupon, eq.([5|) simply asks to apply the zeta-regularized 
product formula once more (at one prescribed point, (V(a) + A)), in order to output the 
absolute-normalized solution ip(a) (or ip'(a)). Thus, eq.(^) acts as the last instruction 
in a procedure to solve the full differential equation ([!]) — for arbitrary A — through 
eqs.(|45"D (still subject to the latter giving convergent iterations, if an effective algorithm 
is wanted). 

As a further prospect, quite general spectral problems can subsequently be tackled in 
principle. E.g., to find the eigenvalues A of eq.([IJ) on the whole line for a general potential 
of even degree, one may proceed to solve eq.(|l|) as just explained but from both ends of 
the g-axis, then match the two resulting values for (i/; 1 /ijj){a) at some finite location a. 

4.2 A numerical test 

We now validate eq.(^B]) with a calculation of the ground state eigenfunction for the 
homogeneous case V(q) = q 4 . The shifted potential is then V a (q) = g 4 +4ag 3 +6a 2 g 2 +4a 3 g, 
and ip(a) is the value of its determinant -D^(A) at the point A = a 4 — Eq, with the eigenvalue 
Eq fa 1.06036209 being part of the input here (while it belongs to the output of an exact 
level calculation for V(q) = q 4 itself [|J). 
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The exact quantization conditions for V a (q) are involving L = 6 distinct chains (except 
at a = 0, where L = 3 by parity symmetry). We only tested calculations of i(j(a) (not 
ip'{a)) for a few values of a > 0, and found the above iteration scheme "C" to converge 
indeed for a < 1.7 (with a contraction ratio per cycle < 0.67). Instabilities made us unable 
to pin down convergent iterations for a > 1.7 (similar comments apply as for scheme "B" 
when t> 2 > +5). Our output points are plotted on fig.3, against the curve produced by a 
standard integration routine and upon which only the global normalization was fitted (a 
posteriori): the results show an overall 4-5-digit agreement. 



4.3 Concluding remarks 

We have reduced the resolution of polynomial Schrodinger equations (|T|) to that of a 
discrete system of selfconsistent exact quantization conditions, eqs. (f45"D , having as un- 
knowns (N + 2) countable sequences of points subject to asymptotic boundary conditions 
consisting of standard Bohr-Sommerfeld formulae (0). Eqs.(^) are supplying explicit 
equilibrium conditions for those semi-infinite and asymptotically tethered chains. 

We also have growing numerical evidence that the so reduced problem is effectively 
solvable in some regions of parameter space through iterations which seem to converge 
geometrically: it thus appears to be a self-stabilizing system, in sharp contrast with the 
original Schrodinger dynamics. This strengthens our following conjecture: (in such cases) 
the relevant equilibrium solution is realized as a fixed point of a contractive mapping, 
which also admits robust finite-dimensional approximations. This results in an overall 
indirect constructive mechanism, where exact quantization formulae explicitly specify the 
mapping only. 



On the darker side, since eqs. (f4~5"|) are very tied to the values of the degree N and 
symmetry order L, they may be ill-suited to transitional regions where one of these num- 
bers jumps (e.g, v 2 — > or + oo in V(q) = q 4 + t>2<? 2 ). Moreover, our tests still span a 
limited range, mostly quartic potentials close to g 4 , a few even sextic ones, and (previ- 
ously |§), homogeneous potentials for their spectra only (but up to quite high degrees). 
Effectiveness for higher N, and extensions to arbitrary complex q and v, to more general 
differential systems, etc., are all conceivable but they remain open issues. 

Finally we argue that, while the integration of eq.([I|) may still be a remote goal using 
quadratures alone, another valuable question is how much the set of admissible integration 
methods needs to be enlarged to reach that same purpose. Our findings provide clues to 
the latter issue: they strongly suggest that zeta-regularized infinite products (of order 
< 1) plus the solving of one type of fixed-point equations (which seem to have contractive 
and other nice properties) are pertinent additional integration tools, which are possibly 
sufficient to the task (already in several cases, at least). 

Acknowledgment: we are grateful to R. Guida (Saclay) for letting us use his com- 
puter program to calculate hundreds of anharmonic oscillator levels with high accuracy. 
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Figure captions 



Fig. 1. Interacting chains for even quartic potentials V(q) = q 4 + v 2 q 2 shown for 
odd spectral sector (i.e., odd quantum number k throughout). The chains occupy their 
equilibrium positions {j e E^ } (points marked with the corresponding t- value). The homo- 
geneous case v 2 = (at center) has full (ternary plus complex-conjugation) symmetry. As 
v 2 — > — oo, the complex chains tend to shadow the resonance spectrum {±{2k + l)iy/— v 2 } 
of the potential — \v 2 \q 2 ; as v 2 — > +oo, they tend to shadow the sequence { — (2k + l)y/v£} 
(£ — 1 chain for the potential v 2 q 2 ). 

Fig. 2. Linear contraction ratios for iteration schemes associated with potentials 
V{q) = q 4 + v 2 q 2 (+: even spectrum; o: odd spectrum). Left: numerical estimates 
for iteration scheme "A" with immediate updating (see main text); although this remark 
is inconclusive, we add that the numerical values of the ratios remain strictly below unity 
(no error estimate at all is implied in our data, but 2-3 digit stability is typically seen). 
Right: numerically estimated moduli of the ratios for iteration scheme "B" which fur- 
ther decouples conjugate-pair interactions (here, various parts of data may correspond to 
different eigenvalue branches, some being negative or in complex pairs). 

Fig. 3. Calculations of the ground-state eigenfunction for the homogeneous quartic 
potential q 4 . +: absolute-normalized data points obtained by iterative exact quantization 
(scheme "C"); in superposition, o: numerical estimates for the corresponding contraction 
ratios. Curve: computer integration of the Schrodinger equation by the NAG routine 
D02KEF; its rescaling factor was the only number fixed by a fit. 
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